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Semel and Skumanich 1^1 proposed a method to obtain the absolute electric current 
density, \Jz\, without disambiguation of 180° in the transverse field directions. The 
advantage of the method is that the uncertainty in the determination of the ambiguity 
in the magnetic azimuth is removed. IHere, we investigate the limits of the calcula- 
tion when applied to a numerical lAHD model'^ We found that the combination of 
changes in the magnetic azimuth with vanishing horizontal field component leads 
to errors, where electric current densities are often strong. Where errors occur, the 
calculation gives \ Jz\ too small by factors typically 1.2 ^ 2.0. 

Numerical MHD Model, Magnetic Field, Electric Current Density 

1 Introduction 

The electric current density is an important measure of the non-potentiality of the magnetic field 
in solar active regions. It has the potential to illuminate the active region eruption process. Vertical 
electric current densities can be calculated from vector magnetic fields measured in the photo- 
sphere. They have been repeatedly obtained ever since photospheric vector magnetic fields first 
became routinely available in the I980s['''5>*^'^'^''''i°'"'i2,i3,i4,i5,i6,i7,i8] ^jj ^Yie&e calculations in- 
volve the resolution of the 180° ambiguity in the transverse field directions, which problem is 
caused by the diagnosis of the magnetic field via the Zeeman effect. Semel & Skumanich (1998)'^' 
developed a method to calculate the absolute vertical electric current density from observed vector 
magnetic fields. Their formulation is especially interesting in that it does not require the 180° dis- 
ambiguation of the transverse field directions and therefore removes one uncertainty in the current 
calculation. 
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The Semel & Skumanich method can be outhned as follows. In SI units, Ampere's law reads 
HqJ = V X B, where J is the current density and /Xq is the permittivity of the vacuum. Since 
/Xo = 47r X 10~ [T m A~ ], then 47r = -^^ — where the magnetic field is measured in G 
and the current density is measured in mA vcT^, and the distances are measured in m. Multiplying 
both sides of the equation by S^, By, and B^By, ByB^, an expression for [(V x B)^]^ is 
obtained 

[(V X B),]2 = Blgl + By^ - Blg^gy sin 2(t> (1) 

which is equivalent to 

|V X B|^ = B^gy - Byg^ (2) 

Here, g^ , gy are given by 

^^"^ y dy 2 " dy 2 ^ dx 

Blgy = Bl \: - -B^By^^ yl- Bl—^ (4) 

ox 2 ox 2 Oy 

where B± = ^JB^. + B"^, the strength of the horizontal magnetic component; B^ = B± cos (p 
and By = B± sin (f) are two perpendicular horizontal components and is the magnetic azimuth. 
All are observable quantities, but [(V x B)^]^ does not vary with either (f)ov (f)+ 180°. 

Recently '^^1, we used Equation (2) to calculate the absolute vertical current densities in two 
morphologically different active regions: a simple sunspot region free of flare activity (NO A A 
10001 on 20 June 2002); and an active region of the quadrupolar configuration with multiple flares 
(including a white fight X3 flarel^^l) and coronal mass ejections (NOAA 10030 on 15 My 2002). 
We found that | \ correlated with vertical magnetic field, B\\,in the form | | = a+6| -By | for both 
regions on large scales. The relationship between | J^j and \B\\ \ needs to be examined with more 
samples of active regions, and observational data from different instnunents. Nevertheless, we 
reaUze that the use of | | has the potential to characterize active regions in their flare productivity. 
The absolute vertical current density is particularly useful for studying the mechanical forces due 
to currents induced in moving material '^^1. 

Equation (2) is limited by the lack of disambiguation of the 180° in the magnetic azimuth. Semel 
& Skumanich^^l pointed out that the calculation "blindly" follows the magnetic field continuity, 
leading to failure when (1) the magnetic azimuth suddenly changes direction; (2) the horizontal 
component simply vanishes paralyzing the continuous magnetic field assumption. In this short 
paper, we explore the limits of the | | calculation with reference to a MHD numerical simulation 
developed by Fan & Gibson f-^ "^'. We will describe the MHD numerical simulation in the next 
section. In section 3, we will present errors of Equation (2), and show statistics over pixels. A 
discussion is given in section 4, and a summary in section 5. 



2 MHD Model 

A 3-dimensional numerical simulation of the coronal response to a rising flux tube was presented by 
Fan & Gibson'^'^1 . The system reaches instability when the flux tube contains a twist of 1 .875 x 27r 
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Fi gure l! The MHD simulations aX t = 30,45,56 in the xy-plane above tlie lower boundary, z = 0.00625-L. The 
background images are the vertical magnetic field, B^. The horizontal magnetic fields are plotted in aiTows: white an'ows 
correspond to Bz < 0, and black arrows correspond to Bz 0. Directions of arrows represent Bx directions, and the 
lengths of the arrows represent Bx strengths. The absolute vertical current densities, | Jz|, are plotted in contours, \Jz\ = 
0.01, 0.04, 0.08, 0.12, 0.15, 0.18, 0.21, 0.23. 
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turns about the tube axis between two footpoints when fully emerged. In this work, we use the 
simulation at three times corresponding to different phases of the flux tube evolution. At t = 30 
(all times are given in arbitrary units), the flux tube axis rises above the lower boundary at a sub- 
Alfvenic speed. At t = 45, the flux tube contains 1.5 x 2it turns between the footpoints above 
the lower boundary. The tube undergoes a significant acceleration with both writhing and rising 
motions, but it stays in a 2-dimensional plane. At t = 56, the flux tube contains 1.875 x 27r 
twists between the two footpoints. The simulation reached the critical point at the onset of the kink 
instabiUty. The complexity of the model is represented by the turns of the flux tube, which grows 
with increasing time. 

The model is built in a Cartesian domain in a box with resolution, 240 x 160 x 200 in the units 
X = [-0.75L,0.75L],y = [-0.5L, 0.5L], z = [0, 1.25L], where L is tiie length of tiie box 
edges. In our work, the "photospheric magnetic field" is the field at ^ = 0.003125L just above 
the lower boundary. Figure 1 shows the 2-dimensional magnetic fields in the model at the three 
times. The vertical magnetic fields, B^, are the background images. The horizontal fields, Bj^, are 
plotted with short arrow bars: white arrow bars correspond to < 0, and black ones correspond 
to Bz ^0. I I are plotted as contours over the magnetic fields. 

3 Errors in the | J^l Calculation 

Errors from the application of Equation (2) are evaluated by comparing the true (| Jz|) with the 
derived current density (| Jz|ss) after excluding the current density on the edges of the simulation 
box (it is apparent that the current densities cannot be calculated on edges of the MHD simulation 
box where Bx{x, y) and By{x, y) are discontinued). | J^l is calculated from 47r = -q^ — 
taking the absolute value. \Jz\ss is calculated with Equation (2), where "ss" represents authors of 
the method, Semel & Skumanich. In this work, the derivatives are replaced with the differences 
between neighboring pixels on the x-y plane: 
dBy _ dB^ _ By{x + Ax, y) - By{x - Ax, y) _ B^{x, y + Ay) - B^{x, y - Ay) 
dx dy ~ 2Ax 2Ay ^ ^ 

where Ax and Ay are the grid unit on the x- and y-axes, respectively. The approximation of 
derivatives to differences gives rise to numerical errors. These errors are introduced to both true 
I J^l and derived current densities | J^l,,.,; but they are squared in the \Jz\ss because the differences 
are carried out for B^, By^B^By between neighboring pixels in Equation (2). The smaller are the 
Ax and Ay, the better are the approximations. 



Figure 2 shows maps of A = | | — | | ss- These maps highlight pixels where large disagree- 
ments between | Jj.| and \Jz\ss occur. Figure 3 shows maps of A J^/ representing the fractional 
errors in derived current density. SmaU values of | A J^, | over weak current densities are given equal 
erroneous impressions to those large values of | A | over strong current densities. 



Table 1 presents quantitative evaluations of Equation (2) when the current densities are calculated 
with models at the times t = 30, 45, 56. The analyses are conducted over pixels selected by two 



4 



Li & Fan Sci China Ser G-Phys Mech Astron \ March 2009 | vol. xx | no. xx | xxxx-xxxx 



t=30 




\\\\\\\ W \ \ \ \ \ W 1 H V W 1 1 I 1 1 1 I J . 



/ / / ^ ^ 
/ / / ^ ^ ^. 
/ / / ^ - - 
/ / / ^ - - 
///-■•- 



- - - - ^ V \ "i, \ 

■ > ^ ^ % N \ WS 
■ ■ ■ . < ^ \ \ \ \ |B 



vTTTTTTTTTTTT 

.r ^ / J J J J I t i I 1 I 



•< ^ ^ < V " ' ■ 




' ' ' ' / f / r ; ; r / / 
t ! t J t 1 t ^ 1 t 7 I t f ^ ^ ^ 

1 I I 1 t m 1 1 1 1 1 1 / r 



Figure 2: Demonstration of Type 1 error. Left: AJ^ = \ Jz\ - \ Jz\ss at three times T = 30, 45, 56 from the top to the 
bottom; The black- and white-most brightness represent [-0.005,0.005]. Small yellow boxes highlight the area where \ Jz \ and 
I Jziss disagree with large uncertainties. Right: Vector magnetic fields within corresponding highlighted yellow boxes, which 
bars and the background images have the same meanings as those in the Fig. 1. The contours represent AJz = \Jz \ — \Jz\bs = 
-0.1, -0.05, -0.005, -0.002, 0.002, 0.005, 0.05, 0.1. Red contours represent AJ^ > and blue contours represent AJ^ < 
0. 
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Figure 3; Demonstration of Type 2 error. Left: AJz/Jz at three times T = 30, 45, 56 from the top to the bottom; 
The black- and white-most brightness represent [-0.5,0.5]. Small yellow boxes highlight the area where large ratios AJ^/Jz 
occur. Right: Vector magnetic fields within corresponding highlighted yellow boxes: background image is the vertical mag- 
netic components; the short bars represent the horizontal magnetic fields. The contours represent AJz = \Jz\ — \ Jz\ss = 
— 0.75, —0.5, —0.2, 0.2, 0.5, 0.75. Red contours represent AJ^/Jz > and blue contours represent AJz/Jz < 0. 
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criteria: |AJj,| ^ 0.005 (Case 1) and lAJ^/J^I ^ 0.5 (Case 2). The percentages of problematic 
pixels in either case are hsted in colunms 1st and 5th of Table 1. The average errors on | J^lss 
calculation are calculated by taking median values of |AJj.| over selected pixels for cases 1 and 
2 (2nd and 6th columns, respectively). The medians of [| J^j, j Jxj.,.,] are presented in the 3rd, 4th, 
and 7th, 8th columns, respectively, over respective sets of pixels defined by either | AJ^j ^ 0.005 
or IAJ2/J2I ^ 0.5. Figure 4 shows the correlation between \ Jz\ (x-axis) and \ Jz\ss (y-axis). All 
pixels are represented by a dot. Pixels where | A | ^ 0.005 are marked by "0" symbols as well. 
Table 1 Statistical Analyses over Pixels 

Case 1: Pixels (| AJ^I ^ 0.005) Case 2: Pixels (|AJ^/ ^ 0.5) 



median median 



Time 


Pixel % 


|AJ,| 




\J-A 


1 -^2 1 SS 


Pixel % 


|AJ,| 




1 t/z 1 SS 


30 


0.308 


8.97 X 10" 


-3 


0.092 


0.049 


3.941 


5.54 X 10-^ 


5.09 X 10-5 


6.13 X 10-5 


45 


0.425 


8.35 X 10" 


-3 


0.193 


0.166 


6.515 


8.42 X 10-5 


5.70 X 10-5 


9.42 X 10-5 


56 


0.675 


7.97 X 10" 


-3 


0.060 


0.047 


7.143 


1.22 X 10-* 


8.05 X 10-5 


1.22 X 10-* 


column 


1st 


2nd 




3rd 


4th 


5th 


6th 


7th 


8th 



4 Discussion 

Within entire simulation planes, we found that the average \Jz\ss and \Jz\ are very close to each 
other. The r.m.s.{\Jz\) are 0.02192, 0.0251 and 0.0211 for models at t=30, 45 and 56, respectively. 
In comparison, the r.m.s.{\Jz\ss) are 0.0168, 0.0244, and 0.0207. They are slightly smaller than 
those I J2 1 by a factor typically 1.01 ~ 1 .03. On the other hand, large differences between | \ and 
\Jz\ss occur on some special pixels which are demonstrated in Cases 1 and 2. 

Numbers of pixels in Case 1 are only a few per thousand pixels (1st column in Table 1). Over 
these pixels, the median | J^l (3rd column) are stronger than average \Jz\ of the entire simulation 
planes by factors 12, 24, and 7 suggesting strong current densities concentrate on these pixels. Ex- 
amination of Table 1 and Figs. 1 and 2 shows that the combination of vanishing horizontal field 
strength and dramatic changes in the magnetic azimuth results in large errors. Examples of prob- 
lematic pixels are marked by yellow rectangles in Fig. 2. On the right hand side, the rectangles 
are enlarged to demonstrate the horizontal magnetic field in all pixels. The vertical magnetic com- 
ponents are the background images, and the A J^, are plotted in contours. This confirms what had 
been discussed by Semel & Skumanichl^l that the method fails on pixels where the magnetic field 
lines are discontinuous. We call errors originated from discontinuous field lines Type 1 error. 

Comparing the 3rd with the 4th columns in Table 1 , the Type 1 error results in that \Jz\ss under- 
estimates the true | by factors 1.9, 1.2, and 1.3. This is also evident in Fig. 4 where the majority 
of "0" symbols fall below the \Jz\ = \Jz\ss line. This can be understood as a result of the "blind" 
calculation of | J^jsg. The dramatic change in the magnetic azimuth causes increasing magnitude 
in \Jz\, but Equations (1) and (2) assume the most close magnetic azimuths between neighboring 
pixels, resulting in smaller | Jzl^s- 

Numbers of pixels in Case 2 are a few percent of the total number of pixels (5th column in 
Table 1). Unlike magnetic fields on pixels in Case 1, the horizontal magnetic fields are neither 
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Figure 4: | J;^ | versus \Jz\s3 at three times, t = 30, 45, 56. The dots represent all pixels within the x-y plane. The "0" 
symbols represent pixels having \ AJz \ > 0.005, most of them represent Type 1 errors. 
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vanishing nor experiencing dramatic changes in the azimuth on pixels in Case 2. The sample 
pixels are demonstrated on the right hand side of Fig. 3 corresponding to the yellow rectangles on 
the AJj/Jj. maps. The median \J^\ (7th column) are only a few hundredth of the average | Jj| 
indicating that the vertical current densities are generally weak on these pixels. Finite differences 
are squared in the Equation (2) and are probably responsible for this kind of error. We call it Type 
2 error. 

The problematic pixels demonstrated in Case 2 are 10 times more numerous than those pixels 
demonstrated in Case 1 (1st and 5th colunms). But the median differences between | J^l and \ Jz\ss 
in Case 2 are more than 100 times as small as those in Case 1 (2nd and 6th columns). These mean 
that the dominant errors of the Equation (2) is Type 1 error. In addition, the median | J^, | in Case 
1 are much stronger than those median | in Case 2 by factors 1807, 3385 and 745 (3rd and 7th 
columns). This means that Type 1 error often occurs over strong vertical current densities. Type 2 
error only becomes prominent where current densities are weak. In real observations. Type 2 error 
can be omitted because it is 2 magnitudes smaller than that of the Type 1 error. 

5 Summary 

Equations (1) and (2), proposed by Semel & Skumanich^^l, provide one more tool to estimate 
the vertical current density from vector magnetic field observations in the observing plane. Their 
method is free of the disambiguation of 180° in the transverse field directions, therefore, removes 
an uncertainty from the calculation. We apply the Semel & Skumanich formula to a MHD 
numerical simulation '^'^l to explore the limits of the method. We identified two types of errors 
to the method. Dramatic changes in the magnetic azimuths with vanishing horizontal fields result 
in Type 1 error. The approximation of the derivatives to differences results in Type 2 error. A 
summary is given below. 

1. Equation (2) accurately measures the vertical current density in a majority of pixels (> 95%) 
in the observing plane. 

2. Type 1 error results from the combination of vanishing horizontal fields and dramatic changes 
in the magnetic azimuth. The errors occur because of an inherit insufficiency in the absolute 
current density calculation proposed by Semel & Skumanich'^l . 

3. Type 2 error results in the approximations of derivatives to differences, which are squared 
in the Semel & Skumanich'^l calculations. Errors become prominent where current densities 
are weak. 

4. Type 1 error is about 100 times as large as Type 2 error. In real observations. Type 2 error 
can be omitted. 

5. When errors occur, the derived current density is smaller than the actual current density by 
factors typically 1.2 ~ 1.9. 

6. The number of unreliable pixels increases with the increasing number of twists in the mag- 
netic flux tube system above the boundary. 
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